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Abstract 

We report the first analytical calculation of the distribution of the 
time headways in some special cases of a particle-hopping model of 
vehicular traffic on idealized single-lane highways and compare with 
the corresponding results of our computer simulation. We also present 
numerical results for the time-headway distribution in more general 
situations in this model. We compare our results with the empirical 
data available in the literature. 
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motion. 
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Significant progress has been made in recent years in understanding the 
dynamical phases of vehicular traffic using particle — hopping models [1]; 
these are analogous to the "microscopic" models of interacting classical par- 
ticles studied in statistical mechanics. The time interval between the arrivals 
at (or departures from) a detector site of two successive vehicles is defined 
as the time headway between this pair of vehicles [2]. The distribution of 
time headways is one of the most important characteristics of vehicular traf- 
fic and the aim of this letter is to calculate this distribution, for a wide range 
of densities of vehicles, in the steady state of the Nagel-Schreckenberg (NS) 
model [3] of traffic on highways. 

In the NS model a lane is represented by a one-dimensional lattice of L 
sites. Each of the lattice sites can be either empty or occupied by at most 
one "vehicle". If periodic boundary condition is imposed, the density c of 
the vehicles is N/L where N(< L) is the total number of vehicles. 

In the NS model [3] the speed V of each vehicle can take one of the 
V max + 1 allowed integer values V — 0, 1, V max . At each discrete time step 
t — > t + 1, the arrangement of N vehicles is updated in parallel according to 
the following "rules". Suppose, V n is the speed of the n-th vehicle at time t. 

Step 1: Acceleration. If, V n < V max , the speed of the n-th vehicle is 
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increased by one, i.e., V n — > V n + 1. 

Step 2: Deceleration (due to other vehicles). If d is the gap in between 
the n-th vehicle and the vehicle in front of it, and if d < V n , the speed of the 
n-th vehicle is reduced to d — 1, i.e., V n — > d — 1. 

Step 3: Randomization. If V n > 0, the speed of the n-th vehicle is 
decreased randomly by unity (i.e., V n — > V n — 1) with probability p (0 < p < 
1); p, the random deceleration probability, is identical for all the vehicles and 
does not change during the updating. 

Step 4 : Vehicle movement. Each vehicle is moved forward so that X n — > 
X n + V n where X n denotes the position of the n-th vehicle at time t. 

We shall present our analytical calculations for V max = 1 and compare 
these with the corresponding results of our computer simulation. As ana- 
lytical calculations are too complicated to carry through for V max > 1, we 
have computed the time-headway distributions for all V max > 1 only through 
computer simulation. 

For the convenience of our analytical calculations, we change the order of 
the steps in the update rules in such a manner that does not influence the 
steady-state properties of the model [4]. We assume the sequence of steps 
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2 — 3 — 4 — 1, instead of 1 — 2 — 3 — 4; the advantage is that there is no 
vehicle with V = immediately after the acceleration step. Consequently, if 
Vmax = 1, we can then use a binary site variable a to describe the state of 
each site; a = represents an empty site and a = 1 represents a site occupied 
by a vehicle whose speed is unity. 

We label the position of the detector by j = 0, the site immediately in 
front of it by j = 1, and so on. The detector clock resets to t — everytime 
a vehicle leaves the detector site. We begin our analytical calculations for 
V max = 1 by assuming that the probability of a time headway t, V(t), between 
a "leading" vehicle (LV) and the "following" vehicle (FV) is given by 

V(t) = t £P(t 1 )Q(t-t 1 ) (1) 

ti=i 

where P(ti) is the probability that there is a time interval t\ between the 
departure of the LV and the arrival of the FV at the detector site and Q(t 2 ) 
is the probability that the FV halts at the detector site for t 2 time steps. 

To calculate P(t\), we have to consider spatial configurations at t — 
for which the FV can reach the detector site within t\ steps. This means 
considering configurations upto the maximum seperation of t\ sites from the 
detector site so that even the farthest vehicle can reach the detector site with 
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t\ hops. The configurations of interest are thus of the form 

(100^__0 10) 71 = 1,2,-..*!. 

n 

The underlined zero implies that we have to find the conditional steady- 
state probability for the given configuration subject to the condition that 
the underlined site (detector site) is empty. This probability II (ra), in the 
2-cluster approximation is given by [4,5] 

n(n)=C(l|0){C(0|0)r" 1 , (2) 

where, C gives the 2-cluster steady-state configurational probability for the 
argument configuration and the underlined imply the conditional, as usual. 
The expressions for the various Cs are given by [4,5] 

C(OjO) = C(0|0) = 1 - V - (3) 

C(0|1) = C(l|0) = f (4) 
C(1\0) = C(0\l) = l (5) 

C(l\l) = C(l|l) = l-- (6) 

c 



where 



V = Y q i 1 ~ V 1 " ^ cd ) » ( 7 ) 



q = 1 — p and d — 1 — c. 

For all configurations with t\ > n, t\ — 1 time steps elapse in crossing n — 1 
bonds (as the last bond is crossed certainly at the last time step). Thus 

P(t 1 ) = J2U(n)q n p t ^ n (8) 

71=1 

Using the expression (2) for IT(n) in equation (8) we get 

P(t l )=C(l\0)q[C(0\Oh + P] t1 ' 1 (9) 
Next we calculate Qfa) by considering all configurations of the form 

(1| 10) 

m— 1 

at t — t\. Here, the underlined part implies that the detector site is occupied 
by the FV and we consider configurations with a queue of (m — 1) vehicles 
ahead of it and the foremost vehicle of this queue has an empty site ahead. 
The steady-state probability of this configuration is given by 

n'(m) =C(1| 11^10) = {C(l|l)} m ~ 1 C(l|0) (10) 

TO— 1 

in 2-cluster approximation. To find all possible configurations of this type 
that contribute to the waiting time of t 2 time-steps, we have to consider 
queue sizes upto t 2 (m < t 2 ). Thus, using expression (10) for II' (m) in 

Q(t 2 ) = J2 Il'(m)q m p t2 - m * 2 - 1 C m _ 1 

TO=1 



we get 



Q(t 2 )=c(i\o)q[qC(i\i)+p} 



«2-l 



(11) 



The time-gap distribution, V(t), can now be calculated using equations 
(1), (9) and (11) and is given by 



where the final result (eq.(12)) was obtained using equations (3)-(6) and y is 
given by eq. (7). 

In the approximate calculation of Q{t?) we ignored the fact that, for 
V max = 1, the LV is certainly present at j — 1 at t — 0. If the LV is 
still at site j — 1 at t > t 1: when the FV is at j = 0, it hinders the forward 
movement of the FV; therefore, we now recalculate V{t) incorporating leading 
corrections. The LV itself may be hindered by its predecessor halting at site 
j = 2 and we get a hierachy of terms. We assume that the probability of 
a m-sized queue ahead of, and including the LV, is given by an expression 
identical to equation (10). Using arguments similar to the one used earlier for 
deriving equation (11), one can show that the probability of the LV hopping 
out of site j — 1 at t — t% is given by Q(t 3 ) (equation (11)). The expression 



V(t) = g 2 C(l|0)C(l|0) 



[C(0\Q)q + pf- 1 - [qC(l\l) + P] 
[C(0\0)q + p}-[qC(l\l)+ P ] 




(12) 
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(9) for P(ti) remains unaffected but the equation (1) gets modified by this 
correction. The expression for the total probability (i.e. V(t)) is now given 
by 

V{t) = ]T {qpt-^PitjSih) +qP(t 1 )R(t 1 ,t-t 1 )} (13) 
ti=i 

where 

S(h) = Q(l) + Q{2) + ■■■ + Q(h) = J2 Q(t 3 ) (14) 

*3=1 

and 

R(h, h) = Q(t 1 +i)p t *- 2 +Q(t 1 +2)p t *~ 3 +- ■ .+Q( tl +t 2 -i) = 2 Q(ts)p tl+ ' 

*3=*1 + 1 

(15) 

The first series on the RHS of equation (13) accounts for the situation where 
the LV hops out of the site j — 1 before the FV arrives at the site j — and 
the second series accounts for the situation where the FV arrives at j = 
before the LV hops out of the site j = 1. The general term in equation (15) 
means that the LV halts for t 3 time-steps (t 3 > ti), the FV being blocked at 
site j = for (t 3 — ti) time-steps, contributing only a factor of Q(t 3 ). After 
the LV hops out of site j — 1 (at t — t 3 ), the FV continues to halt at site 
j = upto t = ti + 1 2 by braking, picking up a factor p*i+*2-*3-i_ 
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Using the expression (11) for Qfa) in (14) and (15) we get 



S(t 1 ) = l-[C(l\l)q+p] t 



(16) 



and 



R(ti,t2) = ^ ) [C(i\i)Q + p] 



pta-l _ [C(l\T)q+p] 



*2-l 



(17) 



So, finally, using (16) and (17) in (13) the total probability distribution 
V(t) can be written as 



V(t) 



+ 



cm 



cm 



cm 

q 2 A C(1|0) 
p-ABC{\\\) "C(l|l 



- 9: 



cmi) 

C(1|0) 



t-i 



(AB) 



t-i 



q 2 A C(1\0) C(1\0) 



p-ASC(l|l) C(0|0) 



p 



(18) 



where A = 1 — gC(l|0), B = 1 — gC(l|0) and equation (18) was obtained 
using equations (3)-(6) and y is given by eq. (7). 

The equation (18) is in excellent agreement with the corresponding results 
of our computer simulation (see figs.l (a) and (b)). Moreover, the distribution 
for the densities c and 1 — c are identical; this is a conseqence of the particle- 
hole symmetry in the problem when V max is unity 

One of the simplest dynamical models of interacting particle systems 
is the so-called asymmetric simple exclusion process (ASEP) [6] which is 
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sometimes regarded as a caricature of vehicular traffic; in this model one 
particle is picked at random and moved forward by one lattice site if the new 
site is empty. In this random sequential update, it is the random picking 
that introduces stochasticity (noise) into the model. Although speed of the 
particles do not explicitly enter into the rules, the effective speed of a particle 
in the ASEP can take only two values, namely, and V max = 1- Since 
mean-field theory is known to give exact results for ASEP, the time-headway 
distribution for ASEP can be obtained from eq.(18) by using the mean-field 
approximations C(l 1 1) = C(0|1) = c = C(l|l) =C(1|0) andC(l|0) = C(0|0) = 
l-c = C(0|l) =C(0|0). 

We present our numerical data for the time-headway distribution in the 
NS model with V max = 5 in fig. 2 as several earlier works have demonstrated 
that this particular choice of V max leads to quite realistic qualitative descrip- 
tions of some other features of vehicular traffic on highways. The qualitative 
features of the distribution in fig. 2 are similar to those in fig.l, except that 
V(t = 0) = but V(t = 1) need not vanish when V max = 5 in contrast to 
the fact that Vit < 1) = 0, irrespective of the density of the vehicles, when 

V m a X 1- 

We have plotted the most probable time-headway T mp as a function of 
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the density in fig. 3; while the curve is symmetric on the two sides of c = 
1/2 when V max = 1 no such symmetry is exhibited for V max > 1. The 
common qualitative trend of variation of T mp is consistent with one's intuitive 
expectation that both at very low and very high densities there are long time 
gaps in between the departures of two successive vehicles from a given site. 

Finally, the time-headway distributions and the trend of their variation 
with density are in qualitative agreement with the coresponding empirical 
data [2]. The time headway distribution in the NS model can be approxi- 
mated well with the forms V(t) oc (t — r) A e~ M * with r = 1 for V max = 1 and 
t = for V max = 5 where the parameters A and /x depend on the density 
of the vehicles [7]. Our results demonstrate that, in spite of being only a 
minimal model, the NS model captures the essential qualitative features of 
the time-headway distribution of vehicular traffic on highways. 
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Figure Captions: 

Fig.l: The time- headway distribution in the steady-state of the NS model 
with V max — 1 (p — 0.5) for the densities (a) c = 0.1(+), c = 0.25(x) and 
c = 0.5 (*). (b) c = 0.9(+), c = 0.75(x) and c = 0.5 (*). The continuous 
curves represent the analytical result while the discrete data points have been 
obtained from computer simulation. 

Fig. 2: The time- headway distribution in the steady-state of the NS model 

with V max = 5 (p = 0.5) for the densities c = 0.10(+), 0.25(x), 0.50(*), 0.75(D), 0.90(b). 

The discrete data points have been obtained from computer simulation while 

the continuous curves are merely guides to the eye. 

Fig. 3: The most probable time-headway plotted as a function of the density 
of the vehicles in the steady-sate of the NS model. The full line corresponds 
to Vmax = 1 and has been obtained from the expression (18); the correspond- 
ing data obtained from our computer simulation have been represented by 
the symbol +. The discrete data points (represented by the symbol x) cor- 
respond to V max = 5 and have been obtained from computer simulation; the 
dotted line joining these data points serves merely as a guide to the eye. 
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